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Abstract 



Static wireless networks are by now quite well understood mathematically through the ran- 
dom geometric graph model. By contrast, there are relatively few rigorous results on the prac- 
tically important case of mobile networks, in which the nodes move over time; moreover, these 
results often make unrealistic assumptions about node mobility such as the ability to make very 
I— ^ large jumps. In this paper we consider a realistic model for mobile wireless networks which 

Cr , we call mobile geometric graphs, and which is a natural extension of the random geometric 

P I graph model. We study two fundamental questions in this model: detection (the time until a 

(-H given "target" point — which may be either fixed or moving — is detected by the network), and 

■y percolation (the time until a given node is able to communicate with the giant component of 

^ the network). For detection, we show that the probability that the detection time exceeds t is 

G exp(— 0(i/ logt)) in two dimensions, and exp(— 0(i)) in three or more dimensions, under reason- 

able assumptions about the motion of the target. For percolation, we show that the probability 
04 that the percolation time exceeds t is exp(— $7(^3+2)) jn all dimensions d > 2. We also give 

^ a sample application of this result by showing that the time required to broadcast a message 

through a mobile network with n nodes above the threshold density for existence of a giant 
component is 0(log ' n) with high probability. 
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1 Introduction 

A principal focus in wireless network research today is on mobile ad hoc networks, in which nodes 
moving in space cooperate to relay packets on behalf of other nodes without any centralized in- 
frastructure. Although the static properties of such networks are by now quite well understood 
mathematically, the additional challenges posed by node mobility have so far received relatively 
little attention from the theory community. In this paper we consider a mathematical model for 
mobile wireless networks, which we call mobile geometric graphs and which is a natural extension of 
the widely studied random geometric graphs model of static networks. We study two fundamental 
problems in this model: the detection problem (time until a fixed or moving target is detected by 
the network), and the percolation problem (time until a given node is able to communicate with 
many other nodes). 

In the random geometric graph (RGG) model [27], nodes are distributed in a region S C M*^ 
according to a Poisson point process of intensity A (i.e., the number of nodes in any subregion 
^ C 5 is Poisson with mean X\A\, where |^| is the volume of A). Two nodes are connected by an 
edge iff their distance is at most r, where the parameter r is the transmission range that specifies 
the distance over which nodes may send and receive information; since the structure of the RGG 
depends only on the product X\Br\ (where Br is the radius-r ball in MJ^) [6], we may fix r so that 
\Br\ = 1 and parameterize the model on A only. We shall take 5" to be a cube of volume n/X (so 
that the expected number of nodes in S is n), and consider the limiting behavior as n — )• oo. 

Clearly, increasing A increases the average degree of the nodes. As is well known, there are two 
critical values of A at which the connectivity properties of the RGG undergo a significant change. 
First there is the percolation threshold A = Ac (a constant that depends on the dimension d), 
so that if A > Ac the network w.h.pj^ has a unique "giant" component containing a constant 
fraction of the nodes, while if A < Ac all components have size O(logn) w.h.p. |27] . Second, at 
the connectivity threshold X = logn, the network becomes connected w.h.p. |19j . The percolation 
threshold A = Ac occurs also in the infinite- volume limit where S" = M'^, in which case the giant 
component is the unique infinite component (or "infinite cluster") with probability 1. These and 
other fundamental properties of RGGs are extensively discussed in the book of Penrose [27j; see 
also ^17j for additional results on thresholds. There are a host of theoretical results on routing and 
other algorithmic questions on static RGGs (see the Related Work section below for a partial list). 
Naturally, most of these consider networks above the connectivity threshold. 

A central feature of many ad hoc networks is the fact that the nodes are moving in space. This 
is the case, for example, in vehicular networks (where sensors are attached to cars, buses or taxis), 
surveillance and disaster recovery applications where mobile sensors are used to survey an area, 
and pocket-switched networks based on mobile communication devices such as cellphones. Such 
networks are also frequently modeled using RGGs, augmented by motion of the nodes. We will 
employ the following model, which we refer to as mobile geometric graphs (MGGs) and which is 
essentially equivalent to the "dynamic boolean model" introduced in [5J in the context of dynamic 
continuum percolation. We begin at time with a (infiniteFl RGG Gq in S" = M'^. Nodes move 
independently in continuous time according to Brownian motion with variance s^; here s is a range 
of motion parameter, which we assume is constant to ensure a realistic model. We observe the 
nodes at discrete time steps (so the displacement of a node in each direction in each time step is 
normally distributed with mean and variance s^). It is not hard to verify that this produces a 
Markovian sequence of RGGs Gq,Gi, . . ., all with the same value of A. Note that, while each Gi 



^We shall take the phrase "w.h.p." ("with high probability") to mean "with probability tending to 1 as n — >■ oo." 

^Passing to infinite volume is a standard device that eliminates boundary effects in a finite region; with a little 

more technical effort the model and results can be extended to finite regions with a suitable convention — such as 



reflection or wraparound — to handle motion of nodes at the boundaries. See, e.g., Corollary 1.3 below. 



is a RGG, there are correlations over time; it is this feature that makes mobihty challenging to 
analyze. 

Once mobility is injected, the questions of interest naturally change from those in the static 
case. For example, connectivity no longer plays such a central role because mobility may allow 
nodes u, v to exchange messages even in the absence of a path between them at any given time: 
namely, u can route its message to v along a time-dependent path, opportunistically using other 
nodes to relay the message towards v. Networks of this kind are often termed "delay tolerant 
networks" [13]. This allows us to focus not on the rather artificial connectivity regime (where A 
grows with n), but instead on the case where A (and hence the average degree) is constant. Keeping 
A constant is obviously highly desirable as it makes the model much more realistic and scalable. 

There are rather few rigorous results on wireless networks with mobile nodes, and those that do 
exist typically either make unrealistic assumptions about node mobility (such as unbounded range 
of motion |18| [TTl |8] or no change in direction [25]), or work in the connectivity regime which, as 
we have seen, requires unbounded density or transmission range [I0l[l3]. (See the Related Work 
section for more details.) In this paper we study two fundamental questions for mobile networks 
assuming only constant average degree and bounded range of motion (i.e., constant values of the 
parameters A and s). 

Results 

Detection. A central issue in surveillance and remote sensing applications is the ability of the 
network to detect a "target" u £ S (which may be either fixed or moving), in the sense that there 
is a node within distance at most r of u. It is well known [26] that, for a static RGG, a fixed target 
can be detected only with constant probability unless the average degree grows with n. In the 
mobile case, we may hope to achieve detection over time with constant average degree (A = 0(1), 
even below the percolation threshold). In this scenario, the detection time, T^et, is formulated as 
the number of steps until a target initially at the origin is detected by the MGG. Recent work of 
Liu et al. [25] shows that the detection time in two dimensions is exponentially distributed when 
the nodes of the network move in fixed directions. In the more realistic MGG model, we are able 
to prove the following result which holds in all dimensions (see Section^: 

Theorem 1.1. In the MGG model with any fixed A and range of motion s > 0, the detection time for 
a fixed target or a target moving under Brownian motion satisfies Pr [T^et ^ t] = exp (— 0(t/logt)) 
for d = 2, and Pr [Tdct > t] = exp (-e(t)) for d>3. 

The constants in the here depend only on A, s and the dimension d. Thus the tail of the 
detection time is exponential in three and higher dimensions, and exponential with a logarithmic 
correction in two dimensions. We note that, as is evident from the proof, this dichotomy between 
two and three dimensions reflects the difference between recurrent and transient random walks 
in I? and Z^ respectively. We also note that the upper bound in Theorem |l.l| holds for arbitrary 
motion of the target (provided it is independent of the motion of the nodes); and the lower bound 
holds for any "sufficiently random" motion of the target. 

We should point out that, for the special case of a fixed target, a slightly stronger version of 
Theorem 1 1.1[ with a tight constant in the exponent, follows from classical results on the "Wiener 
sausage" in continuum percolation. (This was pointed out to us by Yuval Peres [29]; see Related 
Work for details.) However, it is not clear how to extend this approach to the case of a moving 
target. Our proof is elementary and based on an application of the mass transport principle. 

Percolation. A fundamental question in mobile networks is whether a node can efficiently com- 
municate with other nodes even when the network is not connected at any given time. In the MGG 
model, this question may naturally be formulated by considering a constant intensity A > Ac (i.e.. 



above the percolation threshold) and asking how long it takes until a node initially at the origin 
belongs to the giant component (or the infinite component in the limit n — )• oo). We call this the 
percolation time Tperc- It should be clear that the percolation time can be used to derive bounds 
on other natural quantities, such as the time for a node to broadcast information to all other nodes 



(see Corollary 1.3 below). As far as we are aware the percolation time has not been investigated 
before, largely because previous work on RGGs has focused on networks above the connectivity 
threshold. However, it appears to be a fundamental question in the mobile context. 

The detection time clearly provides a lower bound on the percolation time, so we may deduce 
from Theorem |1.1| above that Pr [Tperc > t] is at least exp (—0(t/ log t)) for d = 2 and at least 
exp (— 0(t)) for d > 3. We are able to prove the following stretched exponential upper bound in all 
dimensions d>2 (see Section HI): 

Theorem 1.2. In the MGG model with any fixed A > Ac and range of motion s > 0, the percolation 
time for a node at the origin satisfies Pr [Tperc ^t\= exp (— ri(t ' ' "'"^^)) in all dimensions d> 2. 

Again, the constant in the Jl depends only on A, s and d. There is a gap between this upper 
bound and the lower bound from Theorem |1.1[ We conjecture that the true tail behavior of Tperc 
is exp (-e(t/logt)) for d = 2 and exp {-Q{t)) for d>3. 

Theorem |1.2| is the main technical contribution of the paper; we briefly mention some of the 
ideas used in the proof. The key technical challenge is the dependency of the RGGs Gi over time. 
To overcome this, we partition M into subregions of suitable size and couple the evolution of the 
nodes in each subregion with those of a fresh Poisson point process of slightly smaller intensity 
A' < A which is still larger than the critical value Ac- After a number of steps A that depends 
on the size of the subregion, we are able to arrange that the coupled processes match up almost 
completely. As a result, we can conclude that our original MGG process, observed every A steps, 
contains a sequence of independent Poisson point processes with intensity A" > Ac. (This fact, 
which we believe is of wider applicability, is formally stated in Proposition 4.1 in Section El) This 



independence is sufficient to complete the proof. The slack in the bound comes from the "delay" A. 



To illustrate a sample application of Theorem 1.2, we consider the time taken to broadcast a 
message in a network of finite size. Consider a MGG in a cube of volume n/X (so the expectedj 
number of nodes is n). Since the volume is finite, we need to modify the motion of the nodes to 
take account of boundary effects: following standard practice, we do this by turning the cube into a 
torus (so that nodes "wrap around" when they reach the boundaries) . Suppose a message originates 
at an arbitrary node at time 0, and at each time step t each node that has already received the 
message broadcasts it to all nodes in the same connected component. (Here we are making the 
reasonable assumption that the speed of transmission is much faster than the motion of the nodes, 
so that messages can travel throughout a connected component before it is altered by the motion.) 
Let Tbc denote the time until all nodes have received the message. 

Corollary 1.3. In a MGG on the torus of volume n/X with any fixed A > Ac and range of motion 
■s > 0, the broadcast time Tbc is 0(log ^ ' n) w.h.p. in any dimension d> 2. 

Related work 

There are many theoretical results on routing and other algorithmic questions on (static) RGGs; we 
mention just a few highlights here. The seminal work of Gupta and Kumar [201121] (see also [15j for 
refinements) examined the information-theoretic capacity (or throughput) of such networks above 
the connectivity threshold, i.e., the number of bits per unit time that each node u can transmit 



^The result can be adapted to the case of a fixed number of nodes n using standard "de-Poissonization" argu- 



ments [27j . See the Remark following the proof of Corollary 1.3 in Sectionlsl 



to some (randomly chosen) destination node t„ in steady state, assuming constant size buffers in 
tlie network. Tlie capacity per unit node is 6(n~^/^), wliicfi tends to as n — >• oo, suggesting a 
fundamental limitation on the scalability of such static networks. 

The detection problem has received much attention. In the static case detection is essentially 
equivalent to coverage of the region S, which requires that the network be connected. In the 
absence of coverage, Balister et al. [3] determine the maximum diameter of the uncovered regions, 
while Dousse et al. |12| prove that, for any A > 0, the detection time for a target moving in a fixed 
direction has an exponential tail. (Note that this is not a mobility result as the nodes are fixed.) 

The question of broadcasting within the giant component of a RGG above the percolation 
threshold was recently analyzed by Bradonjic et al. [_7J, who also show that the graph distance 
between any two (sufficiently distant) nodes is at most a constant factor larger than their Euclidean 
distance. Cover times for random walks on (connected) RGGs were investigated by Avin and 
Ercal [T] and Cooper and Frieze j9] , while the effect of physical obstacles that obstruct transmission 
was studied by Frieze et al. [H] . 

The scope of mathematically rigorous work on RGGs with mobility is much more limited. We 
briefly summarize it here. 

Motivated by the fact mentioned above [20] that the capacity of static networks goes to zero 
as n — )• oo, Grossglauser and Tse [18j (see also |11| ) showed how to exploit mobility to achieve 
constant capacity using a two-hop routing scheme. However, these results require the unrealistic 
assumption that nodes move a distance comparable to the diameter of the entire region S at each 
step. El Gamal et al. p!3] study the tradeoff between capacity and delay in a realistic mobility 
model but above the connectivity threshold. Clementi et al. [8] show how to exploit mobility to 
enable broadcast in a RGG sufficiently far above the percolation threshold. However, this result 
again assumes that the range of motion of the nodes is unbounded (i.e., s grows with n). 

As mentioned earlier, the detection problem was addressed by Liu et al. [25], assuming that 
each node moves continuously in a fixed randomly chosen direction; they show that the time it 
takes for the network to detect a target is exponentially distributed with expectation depending 
on the intensity A. Also, for the special case of a stationary target, as observed in [231 |2^ a 
slightly stronger version of Theorem |1.1[ with tight constants in the exponent, can be deduced 
from classical results on continuum percolation: namely, in this case it is shown in [31] that (in 
continuous time) Pr[r(jet > i] = exp(— AKj(t)), where Vs{t) is the expected volume of the "Wiener 
sausage" of length s'^t (essentially the trajectory of a Brownian motion "fattened" by a disk of 
radius r). This volume in turn is known quite precisely |301 13]- 

A model essentially equivalent to MGGs was introduced under the name "dynamic boolean 
model" by Van den Berg et al. [5j, who studied the measure of the set of times at which an 
infinite component exists. Finally, recent work of Diaz et al. [10] in a similar model determines, for 
networks exactly at the connectivity threshold, the expected length of time for which the network 
stays connected (or disconnected) as the nodes move. However, this question makes sense only for 
very large values of A (growing with n) and thus falls outside the scope of our investigations. 

2 Preliminaries 

For any £ > 0, let Bi be the d-dimensional ball centered at the origin with radius i. Similarly, let Qi 
be the cube with side- length £ centered at the origin and with sides parallel to the axes of M*^. For 
any point z G M and set A C M'^, we define z -\- A as the Minkowski sum z + A = {y: y — z ^ A}. 
The volume of a set ^ C M^ is denoted |^|. 

Poisson point processes 

A "point process" is a random collection of points in M'^; for a formal treatment of this topic, the 



reader is referred to [31J. To avoid ambiguity, we refer to the points of a point process as nodes 
and reserve the word points for arbitrary locations in M'^. We are mostly interested in Poisson 
point processes. A Poisson point process with intensity A in a region 5* C M'* is defined by a single 
property: for every bounded Borel set A C S", the number of points in A is a Poisson random 
variable with mean A|^|. We will make use of the following standard properties of Poisson point 
processes: (1) for disjoint sets A, A' C S, the numbers of points in A and in A' are independent; 
(2) conditioned on the number of nodes in A, each such node is located independently and uniformly 
at random in A; (3) [thinning] if each node of a Poisson point process with intensity A is deleted 
with probability p, the result is a Poisson point process with intensity (1 — p)X; (4) [superposition] 
the union of two Poisson point processes in S with intensities Ai and A2 is a Poisson point process 
with intensity Ai + A2. In some of our proofs we will make use of non-homogeneous Poisson point 
processes, whose intensity A(x) is a function of position x G M . In such a process the expected 
number of nodes in a set A is J^ \{x)dx. 

Random geometric graphs 

Fix parameters A, r > 0, and let Sn = Qin/xy/'^ be the cube of volume n/X in R'^. Let H„ be 
a Poisson point process over Sn with intensity A. A random geometric graph (RGG) ^(H„,r) is 
constructed by taking the node set to be the nodes of H„ and creating an edge between every pair 
of nodes whose Euclidean distance is at most r. The parameter r is called the transmission range. 
Since H„ is a Poisson point process, the expected number of nodes in ^(H„,r) is n. 

It is well known [6^,26] that as n — )• 00 the random graph model induced by ^(H„,r) depends 
only on the product A|-Br|. For this reason, we will always fix r = r{d) so that \Br\ = 1 and 
parameterize the model only on A. Note that with this convention, in the limit as n — t- 00, A is also 
the expected degree of any node in ^(H„,r). 

Using a Poisson point process rather than a fixed number of nodes is a standard trick that 
simplifies the mathematics. Most results in this model can be translated to a model with a fixed 
number n of nodes in Sn using a technique known as "de-Poissonization" [27] . 

Many asymptotic properties of ^(H„, r) as n — )• 00 are studied in the monograph by Penrose |27j . 
For example, it is known that A = logn is a threshold for connectivity, in the sense that if A = 
logn + a;(l) then^(H„,r) is connected w.h. p., and if A = logn — a;(l) then^(H„,r) is disconnected 
w.h.p. Another important critical value is the percolation threshold A = Ac (a constant that depends 
on the dimension d); if A > Ac then w.h.p. ^(H„, r) contains a unique "giant" connected component 
with &{n) nodes, while all other components are of size 0(log '^ ~ ' n); on the other hand, if 
A < Ac then w.h.p. all connected components of ^(H„,r) have size O(logn). The value of Ac is not 
known exactly in any dimension d > 2. However, for d = 2 the rigorous bounds 2.18 < Ac < 10.60 
are known [26', Section 3.9], while Balister et al. p] used Monte Carlo methods to deduce that 
A e (4.508, 4.515) with confidence 99.99%. 

Finally, we remark that in the limit as n — )• 00 (that is, when the Poisson point process is 
defined over the whole of M'^) the percolation threshold A = Ac still exists and is characterized by 
the appearance of a unique infinite component (or "infinite cluster") with probability 1 for any 
A > Ac. In this limit the graph is disconnected with probability 1 for any value of A. 

Mobile geometric graphs 

We define our mobile geometric graph (MGG) model by taking a Poisson point process with inten- 
sity A in M at time and letting each node move in continuous time according to an independent 
Brownian motion. We sample the locations of the nodes at discrete time steps i = 1,2, . . . and use 
these locations to define a sequence of random geometric graphs with transmission range r = r{d). 
We base our MGG model on the infinite volume M'^ to avoid having to handle boundary effects on 
the motion of the nodes. Results in this model can be translated to finite regions with a suitable 
convention — such as wraparound or refiection — to handle the motion of nodes at the boundaries. 



More formally, let Ho be a Poisson point process with intensity A over M . We take a parameter 
s > 0, and with each node v £ Hq we associate an independent d-dimensional Brownian motion 
{Wv{i)}i>o that starts at the location of v in Ho and has variance s^ [22]. Now, for any i G Z^, 
we define IIj as the point process obtained by putting a node at Wy{i) for each v G Ho. A MGG 
is then the collection of graphs G = {Gi}i>o where Gi = G{I[i,r) and r = r{d) is fixed so that 
\Br\ = 1. (Note that, as in the static case, fixing the value of r may be done w.l.o.g.) 

It is an easy consequence of the mass transport principle (see below) that each IIj, viewed in 
isolation, is itself a Poisson point process with intensity A. This means that the sequence {nj}j>o is 
stationary and therefore, when viewed in isolation, Gi is a random geometric graph over M*^. Thus, 
for example, if A > Ac then each Gi contains an infinite component with probability 1. 

Mass transport principle 

For two points x, y G M*^ and a time step i > 0, we define fi{x, y) as the probability density function 
for a node located at position x at time to be at position y at time i. Since nodes move according 

II l|2 

to d-dimensional Brownian motion, we have fi{x,y) = j^ — l-)d/2 exp(— 2s'^i ^' 

In some situations, it is useful to regard fi{x,y) as a mass transport function. For example, 
suppose nodes are initially distributed according to a Poisson point process with intensity A in 
a region j4 C M'^; we may view this as a Poisson point process over M*^ with (non- homogeneous) 
intensity (or "mass function") i^o{x) = Alj^-g^}. Using the thinning and superposition properties, 
it is easy to check that the distribution of the nodes at time z is a Poisson point process with 
intensity fj(y) = J^a t^oix) fi{x,y)dx. This interpretation can be used, for example, to show that in 
a MGG IIj is a Poisson point process with intensity A for all i, as claimed above. 

3 Detection time 



In this section we prove Theorem 1.1 We consider the detection time for a node u initially placed 
at the origin independently of the MGG G. We say that a node v G G detects u at time step i if 
the distance between u and v at time step i is at most r, and we define T^et as the first time that 
u is detected by some node of G. Our goal is to derive tight bounds for the tail of Pr [T^ct > t]. 
In the proof we consider the cases where u is either non-mobile or moves according to Brownian 
motion with variance s^. We discuss some extensions at the end of the section. 

It will be convenient to restrict attention to the nodes of G that are initially inside the cube Ql, 
where L is a suitably chosen parameter. We define TfictiQL) as the first time a node initially inside 
Ql detects u. Note that clearly Pr[TAet{QL) >t] > Pr [Tdet >t] = limL^^Pr[Tdet{QL) > t], 
where the limit exists since Pr [T^etiQi] ^ t] is monotone and bounded as a function of L. We let 
X = {xq,xi, . . . ,xj_i) be the locations of u in the first t steps. The following lemma relates the 
tail of Tdct(QL) to the tail of an analogous random variable for a single random node in Ql. 

Lemma 3.1. We have Pr [TdetiQi] > t | X] = exp {-XL'^ Pr [r < t | X]) and Pr [TdctiQL] >t]> 
exp (— AL Pr [r < tfj, where r is the first time that a node initially located u.a.r. in Ql detects u. 

Proof. Let A^ be the number of nodes inside Ql at time 0. Each of these nodes is initially located 
uniformly at random inside Ql, and the motion of each node does not depend on the locations of 
the other nodes. If we fix a given value for X, then the first time that a given node v detects u does 
not depend on the other nodes of G and is distributed according to the conditional distribution 
of r given X . (Note that this is not true if X is not fixed but random, because the random motion 
of u makes the relative displacements of the nodes of G with respect to u dependent.) Therefore, 
conditioning on N and X, we have Pr [T^ct^QL) >t\N = n,X]= Pr [r > i | X]"", which yields 

Pr [TdotiQi] > i I A] = Ejv [Pr [r > t I xf] = exp (-AL'^Pr [t < t \ X] 



where we use the notation E^r [•] to denote expectation with respect to the random variable A^, and 
the last equality holds since A^ is Poisson with mean AJQlI = AL'^. For the lower bound, we appeal 
to Jensen's inequality to obtain 

Pr[Tdet(QL)>t] = Ex[Pr[Td,t{QL)>t\X]] 

> exp (-XL'^ Ex [Pr [t < t \ X]]) 

= exp('-AL'^Pr[r <t]V 

which completes the proof. D 

We now proceed to derive upper and lower bounds for Pr [t < t]. Let w be a node initially 
located u.a.r. in Ql- For time steps ii < 12, let M{ii,i2) be the expected number of time steps 
from ii to Z2 at which v detects u. We bound M(0, t — 1) as follows. 

Lemma 3.2. Let t G Z+ and X = (xo,a;i, . . . ,xt-i) G M be arbitrary. There exists a constant 
c = c{d) and Lq = Lo(t, maxj ||xj||2) such that ct/L'^ < M(0, t — 1) < t/L'^ for all L > Lq. 

Proof. We use the mass transport principle. We assume the initial intensity i^o{z) = Al/^gg^i and 
let fj(z) be the intensity at z S M'^ at time i, i.e., ^{{z) = f^^ '^o{y)fi{y, z)dy for i > 1. At any time 
i, the probability that v detects u is given by the ratio between the amount of mass inside Xi + Br 
and the total amount of mass XL'^. Noting that for i = this ratio is ynr = 1/-^'^) we can write 
M(0, t - 1) as 

M(0,t-1) = 



1 


^ ^ XLd 

i=l 


z)dz 




1 


^ >•!'■' i^JsjQ, 


Mi{y,Xi 


+ z)dydz 


1 


^ j = l J Br Jxi+Z + QL 


y')dy'dz, 



(1) 

where the last step follows from the translation-invariance property of fi . Since \Br\ = 1 , we obtain 
the upper bound from /3,^+^+Q^ fi{0,y')dy' < f^a fi{0,y')dy' = 1. 

For the lower bound, let Ai = B^n^. We assume that L is sufficiently large such that Ai C 
Xi + z + Ql ■ Then replacing the integral over Xi + z + Ql in ([l]) by an integral over Ai gives the 
lower bound 

1 l^/"/" ,. 1 1^ lAiWBJ , . ct 



m^^-^)>-T^^T^Y.jJj^M)dy'dz>j, + l,Y^AM^^ 



^^, J Br J A. - .- -. (2.,2,)./2 



Ld' 



where we used the fact that \\y'\\2 < sy/i + 1 for all y' £ Ai, and the result holds for some constant 
c = c{d). D 

Our goal is to write M{ii,i2) conditioning on r. Let M'(y, 21,22) be the expected number of 
time steps from ii to 12 at which v detects u given that the relative location of v with respect to u 
at time ii — 1 is y. The next lemma gives lower and upper bounds for M'{Yi,i + 1, i + 1). 



Lemma 3.3. Leti G Z^ he arbitrary. There exists an integer t^ = tQ{d, s) such that for allt > to the 
following holds. There exist functions mi{t) andm2{t) such thatmi{t) < M'{Yi,i + l,i + t) < m2{t) 
uniformly over Yi ^ Br. Moreover, there are constants ci = ci{d) and C2 = C2{d) such that 

u\ ^ / "^1 logt/s"', ford = 2 ^ / C2 logt/s'^, for d = 2 

^^(*)^[ c,/s\ ford>?, "^2(t)<| ^^/^.^ ^^^^>3 

The hounds for mi{t) hold both for the case where u does not move and for the case where u moves 
according to Brownian motion with variance s^. The hounds for m2{t) hold uniformly over X. 

Proof. For any j £ [l,t], let Ij be the indicator random variable for the event that v detects u 
at time i + j, assuming that at time i u is located at the origin and v is located at YJ. Clearly, 
M'{Yi, i + l,i + t) = X]j=i E [Ij]. Recall that Xj+j is the location of u at time i + j. Hence, 

E [Ij] = / ,^ o .,,/o exp w-^, dz < 



(27rs2jy/2 ^ V 2s2j ; - (27rs2jy/2 (27rs2jy/2 



The upper bound follows by setting m2{t) = '}2j=i 7^ — 2-\d/2 ■ Note that this upper bound holds for 
arbitrary X. 

Now we derive the lower bound. We use the fact that 1^,2; G B^. If u is non-mobile, then 
||xi+j||2 = \\xi\\2 = (recall that we assume Xi to be the origin) and from the triangle inequality 

we obtain \\xi+j + z-Yi\\2 < \\xi+j\\ + ||2;||2 + ||li||2 < 2r. Thus, E [/,■] > (2^g2"j)d/2 exp (-frjj- We 
take jo to be the smallest integer such that jo > r"^ /s"^, set mi{t) = J2j=jo Jo — 2^-1^/2 exp(— 2), and 
the result follows since t is sufficiently large with respect to jq. 

If u moves according to a Brownian motion with variance §2 ^ ^e average over Xj+j to get 

1 / ||'r-',--l-^ V-l|2\ 1 / ll-r-i-lP 

-L / ll-'^i+J 1^ -'«ll2\ / Il'^«+jll2 



Let a > be a constant and let Aj = -B^syj_2r- We set a > 4 so that as^/J — 2r > as^/j/2 
for all j > jo = [r2/s2]. We integrate over Aj instead of M'^ and then use the simple bounds 
\\xi^j + z — Yi\\2 < as^/J and ||xi+j||2 < as^/j for all Xj+j € Aj and z, Yi G i?r to obtain 

E /,• > ^'' ^' exp (-a2) > — -— TTT, 
^ ^^ - (27rs2j)<^ ^^ ^ - (s2j)'i/2' 

for some constant a' = a'{d). Now, we set mi{t) = Y1,j=jq i 2°w/2 ^'i^d the result follows since t is 
sufficiently large with respect to jo- □ 

Remark: The bounds for M'(Yi, i + 1, i + t) change substantially from d = 2 to d > 3, reflecting the 
dichotomy between recurrent and transient random walks in W? and Z^: v returns to a neighborhood 
of u infinitely often for d = 2 and only finitely often for d> 3. (Note that M' measures the expected 
number of returns of u to a neighborhood around t> in a given time interval.) 

We now use Lemma 3.3 to derive upper and lower bounds for Pr [r < t]. 
Lemma 3.4. Let the functions nii and m2 be as in Lemma\3.3[ For any constant a > we have 



^(0'*-l)<Pr[^<i]<^(0'(l+«)*-l) 



1 + m2{t) l + mi(at) 



Proof. We apply the straightforward equation 

t-i 



M{0,t-l) = Y,Pr[T = i]{l + EY,[M'{Yi,i + l,t-l)]), 



i=0 



where the random variable Yi denotes the relative location of v with respect to u given that t = i. 
Note that Yi G Br since the condition t = i implies that the distance between v and u at time i is 



at most r. Using Lemma 3.3 we obtain 
t-i 



M{0,t - 1) <^Pr [t = i] {I +-Ey, [M'{Yi,i + l,i + t)]) < Pr [r < t] (1 + m2(t)) . (2) 



i=0 



Also, since M(0, t — 1) is non-decreasing with t, we can take an arbitrary constant a > and use 



the fact that (1 + a)t — 1 > i + at for all i G [0, t — 1] together with Lemma 3.3 to write 

t-i 
M(0,(l + a)t-l) > Y^Pr [t = i] {1 + Ey^ [M' {Yi, i + 1, {I + a)t- I)]) 

i=0 

t-i 

> ^Pr [r = f] (1 + Ey, [M'{Yi,i + l,i + at)]) 

j=0 

> Pr[T <t]{l+mi{at)). (3) 

The proof is completed by reorganizing ^ and ([s]) . D 

We are now in position to conclude the proof of Theorem |l.l[ Plugging Lemmas 3.2 and 3.3 into 



Lemma 3.4, and using Lemma 3.1, we obtain a constant to = to{d,s), and constants ci, 02,03,04 



depending only on d, such that for all t > tQ and sufficiently large L, 

exp (- , ''^']' ) < Pr [Tdet(QL) > t] < exp (- , ""f']' ] (4) 

V S"' + C2logty V S^ + C4logt/ 

for d = 2, and 

„d 



exp i^-ciXs^'tj < Pr [TdctiQi) > t] < exp ^-C3XsHj (5) 

for d>3. Theorem 1 1 . 1 1 then follows by taking the limit as L — )• 00. 

Remark: As should be clear from the proof, the upper bounds in Q and (pj) hold for arbitrary 
locations of u as long as u moves independently of the locations of the nodes of G. The lower bounds 
also hold in more generality: e.g., if u moves according to Brownian motion with variance s' 7^ s^, 
or indeed with any motion that has sufficiently large "variance" in all directions. (Specifically, 
the lower bounds hold if the density /j'(-, •) for the motion of u after i steps satisfies the following 
property: there exist positive constants ai = ai{d) and 02 = a2{d) such that fl{x,x + z) > ai/i'^''^ 
for all z £ B^ n, X G M , and i G Z"*".) On the other hand, adding a random drift to the nodes in 
G can change the detection time substantially: if each node v oi G moves according to Brownian 
motion with drift fiy and variance s^, where the fiv are i.i.d. random variables, then our proof can be 
adapted to show that under mild condition^ on the distribution of fiy, Pr [T^ct ^ ^] = 6xp (— 0(t)) 
in all dimensions d>2 for arbitrary locations of u. We omit the details. 



''Note that this statement cannot hold in full generality; if all nodes of G have the same drift, then this is equivalent 
to the case without drift up to translations of V^ . 



4 Percolation time 



In this section we prove Theorem 1.2 We consider a MGG G with density A > Ac (i.e., above 
the percolation threshold), and study the random variable Tperc defined as the first time at which 
a node u initially placed at the origin independently of the nodes of G belongs to the infinite 
component of G. We derive an upper bound for the tail Pr [Tperc > i] as t — )• oo. 

We begin by stating a proposition that will be a key ingredient in our analysis. We consider 
a large cube Qk C W^ and tessellate it into small cubes called "cells." The proposition says that, 
if all cells have sufficiently many nodes at a given time i, then at time i + A for suitably large A 
the point process induced by the location of the nodes contains a fresh Poisson point process with 
only slightly reduced intensity inside a smaller cube Qk' ■ We believe this result is of independent 
interest. With this in mind, we state the proposition below for a slightly more general setting than 
is needed here. Its proof is deferred to the end of the section. 

Proposition 4.1. Fix K > I > and consider the cube Qk tessellated into cells of side-length L 
Let IIq be an arbitrary point process at time that contains at least f3£ nodes in each cell of the 
tessellation for some /? > 0. Let IIa be the point process obtained at time A from Uq by allowing 
the nodes to move according to Brownian motion with variance s^. Fix e E (0, 1) and let 'E be a 
fresh Poisson point process with intensity (1 — e)/3. Then there exists a coupling ofE and IIa o-nd 
constants ci, C2, C3 depending only on d such that, if A> %^ and K' < K — C2Sy/Aloge^^ > 0, the 
nodes ofE are a subset of the nodes of Ha inside the cube Qk' with probability 1 — 73- exp(— cse^/?^'^). 



Now we proceed to the proof of Theorem 1.2 We first take a sufficiently small parameter ^ > 
such that (1 — ^)^A > Ac. (This is always possible as we are assuming A > Ac.) In what follows, we 
omit the dependencies of other parameters on A and ^ as we are considering them to be fixed. 

Let Hi be the event that u does not belong to the infinite component at time i. Then, the 
event {Tperc > t} is equivalent to ni=o-^*- ^^ define an integer parameter A > 1 and consider 
the process obtained by skipping every A time steps. (To simplify the notation we assume w.l.o.g. 
that t/A is an integer.) In other words, instead of looking at the event 0^=0 -^* ^^ consider the 
event ni=o ^Ai, which we henceforth denote by Hf Since the occurrence of the event {Tperc > t} 
implies %t we have Pr [Tperc > i] < Pr \}it\- Our goal in introducing A is to allow nodes to move 
further between consecutive time steps; we will choose the value of A later. 

Let C = C{d) > 1 be a sufficiently large constant and fix L = Ct(l + s). We will confine our 
attention to the cube Q2L- We take a parameter £ > and tessellate Q2L into cubes of side-length l 
(see Figure [Ijl^a) ) . We refer to each such cube as a "cell." Later we will tie together the values of (. 
and A, and will choose £ to optimize our upper bound for Pr {Ht\. For the moment we only assume 
that the tessellation is non-trivial in the sense that both ^ and L/(. are (^(1) as functions of t. 

For each time step i, the expected number of nodes inside a given cell is A^*^. We say that a cell 
is dense at time i if it contains at least (1 — ^)Xi nodes, where ^ is as defined earlier. Let Di be 
the event that all cells are dense at time i, and let Dt = HiLo ^i- "^^^ lemma below shows that 
T>t occurs with high probability. 

Lemma 4.2. With the above notation, Pr [Vt] > 1 - ^ exp [-^'^M'^/i) . 



Proof. At any given step i, by a standard large deviation bound for a Poisson r.v. (cf. Lemma A.l ), 
a cell has more than (1 — £^)Xi'^ nodes with probability at least 1 — exp(— ^^A^'^/2). The proof is 
completed by taking the union bound over all {L/l) cells and t/A time steps. D 

Recall that Xj is the location of u at time i. Define Ei as the event that Xi is located inside 
Q^/3 at time i, and let £t = fliLo -^-Ai- The next lemma bounds the probability that u never 
leaves Ql/s- 
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Figure 1: (a) The cubes Q2L and Ql and the tessellation of Q2L into cubes of side- length i. (b) The 
cube Ql/3, the locations xq and Xi of node u at time steps and i respectively, and the cube Si. 



Lemma 4.3. There exists a constant c = c{d) such that Pr [£t] > 1 — exp {—ct) . 

Proof. We fix a time step i and then apply the union bound over time. The event Ei corresponds 
to u not moving a distance more than L/6 in any dimension. Therefore, 



Pr[Ei] = [ fi{Q,xi)dxi 



Ql/3 



1-2 



L/6 



v2vrs^ 



exp 



2s2? 



dy 



> l-2d 



L/6 



v2vrs^ 



exp 



2sh 



dy. 



Then we use a standard large deviation bound for the Normal distribution (see Lemma A. 2) to 
conclude that 

Pr \Ei\ > 1 y==^ exp 



/27rL ^ V 72s^i 
Since the bound above decreases with i we can conclude that 



Pr [£t] > 1 



t l2ds^/t 



A y/2^L 

and the result follows from A > 1 and L > st. 



exp 



72sH 



D 



For each time step i, we define St to be the cube Ql/z shifted randomly so that Xi (the location 
of u at time i) is uniformly random in Si (see Figure [lib)). A crossing component of Si is a 
connected set of nodes within Si that contains a path connecting every pair of opposite faces of Si. 
(A path connects two faces of Si if each face is within distance r of one of its endpoints.) 

For each z, let Ki be the event that all the crossing components of Si are contained in the infinite 
component at time i. (For definiteness we assume that Ki holds if Si has no crossing component.) 
Let JCt = HiLo -^Aj- The next lemma follows by a result of Penrose and Pisztora [28, Theorem 1]. 

Lemma 4.4. For any A > Ac, there exists a constant c = c{d) such that Pr [/Ct] > 1 — exp(— ct'^"^). 

Proof. By stationarity we know that Pr [Ki] is the same for all i. For any fixed i, |28L Theorem 1] 
gives that Pr [Ki] > 1 — exp(— c'L*^"^) for some constant c'. (In fact, [281 Theorem 1] handles an 
event more restrictive than Ki, which among other things considers unique crossing components 
of Si.) Using the union bound we obtain Pr [/Cj] > 1 — (i/A) exp(— c'L*^"^) and the result follows 
since L > t. D 
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We now proceed to derive a bound on Pr \Ht\- We take H[ to be the event that u does not 

belong to a crossing component of Si at time i, and define %[ = HiLo -^Ai- Note that H[ is a 
decreasing event, in the sense that if H[ occurs then it also occurs after removing any arbitrary 
collection of nodes from the MGG G. Clearly 'Htr\K.t ^T~i[r\K.t- By elementary probability, 

Pr [Ht] < Pr [n[ n Vt I £t] + Pr [V^] + Pr [S^ + Pr [/C?] . (6) 

Note that we use /Ct only to replace Tit by ?^^ in ([6]) ; this helps to control the dependencies among 
time steps, since Ti'^ is an event restricted to the cubes Si while Tit is an event over the whole of M'^. 
We use £t only to ensure that Si C Ql, which allows us to focus on the portion of G inside Q^. 
Note that £t is independent of G, so this conditioning does not affect G. 

Now we set A = [C^£^/s^], where C is the constant in the definition of L. The main step in 
our proof is the lemma below. 

Lemma 4.5. Let A > Ac and ^ > be such that (1 — ^)^A > Ac- Let G be large enough in the 
definition of L and A. There exist constants c = c{d) and to = tQ{d) such that, for all t > to, we 
have 

Pr [n't n Vt I £t] < exp (-ct/A) . 

Proof. We start by writing 

t/A-l 

Pr [n't nvt\£t]< n P''[-^A, I 'H'^i^-l) n v^h-d n St]. (7) 

(Here, for notational convenience, we assume that n'_^ n P-a H <?i = £f) 

We now derive an upper bound for Fr[H'^- \ n'^,^_^^ l^^A(j-i) ^^t\- We start with a high level 
overview of the proof. Let ^^(i-i) be the (not necessarily Poisson) point process obtained from the 
nodes of ^/\(i-i) (the MGG at time A(i — 1)) under the condition n'^u_i\ n ^a(j-i)- Note that 
'^A(i-i) is conditioned only on events that occur between time and time A(i — 1); therefore, the 
motion of the nodes of ^^(i-i) from time A(i — 1) to Ai is independent of the condition. Since all 



cells are assumed dense at time A(z — 1), using Proposition 4.1 we can construct an independent 
Poisson point process S'^(j_i) and couple it with Ha(j-i) so that at time Ai the nodes of S'^j in Ql 
are a subset of the nodes of Ha^. Moreover, we can ensure that H^^ has intensity larger than Ac in 
Ql , and thus conclude that u will belong to a crossing component of S/\i with constant probability. 
Using this, we can upper bound each term of the product in ([7]) by a constant strictly smaller 
than 1, which gives Pr \n't H Vt \ St\ < exp(— ct/A) for some constant c = c{d). 

Turning now to the details, we can invoke Proposition |4.1| with /3 = (1 — ^)A, e = ^, K = 2L, 
and K' = L to obtain that, conditioned on n'^,■_^^ n ^a(j-i); at time Ai the nodes of the MGG 
in Ql contain a fresh Poisson point process with intensity (1 — ^)^A > Ac with probability at least 
1 — exp (— c'^^A^*^), for some constant c'. Since H'^^ is a decreasing evenlrl if we define H'^^^ as the 
event H'^- restricted to the nodes of the fresh Poisson point process, then H'^- C H'^- and 



Pr 



nAi 



Ai 



^Ali-i) n ^A(*-i) n £t] < Pr [n"j^,] + exp (-c'e2^^d\ ^ 



as n"^^ does not depend on the condition. Since the intensity of the fresh Poisson point process 
is (1 — ^)^A > Ac, p?^ Theorem 1] implies that, with probability 1 — exp(— c"L''~^) for some 
constant c", a constant fraction of the volume of S'aj is within distance r of at least one node in 



""Note that we defined H'^^ in terms of crossing components precisely in order to make it a decreasing event; 
otlierwise, we could have defined it in terms of the largest connected component of Sai- 
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a crossing component of Saj- Since at time Ai, u is located uniformly at random inside S/^i, u 
belongs to a crossing component of Sai with probability at least c'" > 0, where c'" is a constant. 
Hence we have 

t/A-l 
Pr [n't nVt\£t]< n [l - c'" + exp(-c"L'^-i) + exp (-c'e^A^*^)" . 

Since L and £ go to infinity with t, for t sufficiently large each factor in the above product can be 



made strictly smaller than 1, which concludes the proof of Lemma 4.5 D 

Finally, we plug Lemmas 4.2-4.5 into Q and obtain the following upper bound on Pr [Ht]'- 

Pr [Tperc > t] < Pr [Ht] < exp {-ct/A) + exp (-ce"^) + exp {-ct) + exp(-ct'^-^). 



(Here c = c{d) is a generic constant.) In order to minimize this upper bound we choose i so that 
i'i = e(t/A), which yields 

Pr [Tperc >t]< exp (-ctd2 



for all sufficiently large t, where c is a constant depending on A, s, and d. This completes the proof 



of Theorem 1.2 It remains to go back and prove Proposition 4.1 



Proof of Proposition 4.1 



We will construct H via three Poisson point processes. We start by defining Hq as a Poisson point 
process over Qk with intensity (1 — e/2)/3. Recall that Hq has at least pi'^ nodes in each cell of 
Qk- Then, in any fixed cell, Hq has fewer nodes than Hq if Hq has less than /3£'^ nodes in that 
cell, which by a standard Chernoff bound (cf. Lemma A.l[) occurs with probability larger than 



1 - exp (- '' (1-^/2)/^^" (1 _ e'/s)") for e' such that (1 + e')(l - e/2) = 1. Since e G (0, 1) we have 

e' G (e/2,1), and the probability above can be bounded below by 1 — exp (— ce^/^^'^) for some 
constant c = c{d). Let {Hq ^ Hq} be the event that Hq has fewer nodes than Hq in every cell of 
Qk- Using the union bound over cells we obtain 

Pr[H'Q ^ Ho] > 1 - 1^ exp(-ce2/3^^). (8) 

If {Hq :< Hq} holds, then we can map each node of Hq to a unique node of Hq in the same 
cell. We will now show that we can couple the motion of the nodes in Hq with the motion of their 
respective pairs in Hq so that the probability that an arbitrary pair is at the same location at time 
A is sufficiently large. 

To describe the coupling, let v' be a node from Hq located at y' S Qk, and let v be the pair of 
v' in Hq. Let y be the location of v in Qk, and note that since v and v' belong to the same cell we 
have \\y — y'\\2 < vdi. We will construct a function g{z) that is smaller than the densities for the 
motions of v and v' to the location y' + z, uniformly for z G M . That is, 

r „ / / / N „ / ; M 1 / niaxjllzlln, llv' + z — w||n|\ , , 
9{z) < mm{/A(y', y' + .), fA{y, y' + z)} = j^^^^j^^, exp (^ ^^^f^ ^^ J (9) 



for ah zeW^. 
We set 



1 / {\\zh + ^dlf \ 
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Note that this definition satisfies ([9]) since by the triangle inequaUty Hy' + z — y||2 < ||y' — ylb + H-^lb 
and \\y' — y\\2 < vdL Define il> = 1 — J^^g{z)dz. Then, with probability 1 — ■0 we can use the 

density function j^ to sample a single location for the position of both v and v' at time A, and 
then set Hg to be the Poisson point process with intensity (1 — 'i/')(l — £/2)/3 obtained by thinning 
Hq (i.e., deleting each node of Hq with probability ip). At this step we have crucially used the fact 



that the function g{z) in (10) is oblivious of the location of v and, consequently, is independent of 
the point process Ho. (If one were to use the maximal coupling suggested by ([9]), then the thinning 
probability would depend on Ho, and Hq would not be a Poisson point process.) 

Let H^ be obtained from Hq after the nodes have moved according to the density function ^^ . 
Thus we are assured that the nodes of the Poisson point process r!'^ are a subset of the nodes of 
IIa and are independent of the nodes of XIq, where Ha is obtained by letting the nodes of Ho move 
from time to time A. 

The next lemma shows that if A and K — K' are large enough, then the integral of g{z) inside 
the ball -B(_ft-_A-')/2 is larger than 1 — e/2. (We are interested in the ball B(^j^_j^iy2 since for all 
z G Qk' we have z + B(^k-K')/2 C Qk-) 

Lemma 4.6. If A > c^^n o-nd K — K' > c'syAloge"^ for large enough c,c' , we may ensure that 
L g{z)dz > 1 - e/2. 

Proof. Since g{z) depends on z only via ||z||2, we integrate over p = \\z\\2 and let ap'^~^ be the 
surface area of Bp, which gives 

9i^)dz = / ,^ 9AN^/9 exp ^-^-r dp. 



B,.-.0/. ^0 {2.s^Ar/^ ^V 2.^A 



Now we change variables to p' = ^ y^ and set b = ^^^ = e\/d/c and K" = - — 7= + 5 to obtain 

^ ^ sVA sVA V / 2sVA 

/ 9{z)dz = -^ r [p' - df-' eM-p'^mdp'. (11) 

Let /i(5) = /^°°(/9' — 5)^^"^ exp(— p' /2)dp' . We will apply the Taylor expansion of /i(5) around 5 = 0. 
Note that , ?j/2 ^(0) = 1) ^"^^ the derivative of h{5) is 

h'{5) = -{d - 1) / {p' - SY'"^ exp{-p'^/2)dp'. 

In particular, the derivative increases with 6, and by Taylor's theorem, h{6) > /i(0) + 6h'{0). 
Therefore, we have 

« .h{6) > -^ {h{0) + 6h'{0)) = 1 - 6-"'''^^^ 



(27^)^/2 ^ ' - (27r)^/2 ^ ^ ' ' ^ " (27r)'^/2' 

Note that 6 "" ^^2 depends on the dimension only and can, for example, be made smaller than e/4 



for sufficiently large c. Using equation (11) we get 



/. 



B{K-K')/2 



j'OO 

g{z)dz > 1 - 6/4 - j-^^ j^^^ {p' - 5)'-' exp(-p'V2) V. (12) 
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Now note that p"^ ^ exp(— p' /2) < c" exp(— p' /3) uniformly for p' G [0, oo), where c" is a constant 
depending only on d. Thus we have 

/,. W^^^' - ^^'" -M-p"/2)dp' < ^-^^ j^^^ eM-p"mdp' < c'" eM-K"'/3), 
for a constant c'" depending only on d. For sufficiently large c' (and thus K"), the right hand side 



in ( 13 ) can be made smaller than e/4. Plugging this into ( 12 ) yields the lemma. D 



When {Hq ^ IIo} holds, H'^ consists of a subset of the nodes of Ha- Note that H'^ is a non- 
homogeneous Poisson point process over Qk- It remains to show that the intensity of S^ is strictly 
larger than (1 — e)/3 in Qx' so that H can be obtained from H^ via thinning; since H'^ is independent 
of Ho, so is H. 

For z E M'^, let p{z) be the intensity of H'^. Since Hq has no node outside Qk, we obtain for 
any z G Qk' , 

p{z) > (1 - i^){l - e/2)/3 /" ^T^^"^ = ^^ " '/^^^ / ^^''^'^''' 



where the inequality follows since z + Btx-K')/2 C Qftr foi' ^^ -^ £ Q.K''- From Lemma 4.6 we have 



Jb g{x)dx > 1 — e/2. We then obtain p{z) > (1 — e/2)^/3 > (1 — e)/3, which is the intensity 

{K~K )/2 

of H. Therefore, when {Hq ^ IIo} holds, which occurs with probability given by ([8]), the nodes of 



S are a subset of the nodes of IIa, which completes the proof of Proposition 4.1 

5 Broadcast time 



In this section we use Theorem 1.2 to prove Corollary 1.3 for a finite mobile network of volume 
n/X. 

We may relate the MGG model on the torus to a model on M as follows. Let Sn denote the cube 
Q(n/A)i/d- The initial distribution of the nodes is a Poisson point process over M'^ with intensity A 
on Sn and zero elsewhere. We allow the nodes to move according to Brownian motion over M as 
usual, and at each time step we project the location of each node onto Sn in the obvious way. 

Now let t = Clog ^ ' n for some sufficiently large constant C = C{d). The proof proceeds 
in three stages. First, we show that for any fixed i G [0,t — 1], the giant component of Gi has at 
least one node in common with the giant component of Gj+i. This means that, once the message 
has reached the giant component, it will reach any node v as soon as v itself belongs to the giant 
component. Thus we can bound the broadcast time by (twice) the time until all nodes of G have 
been in the giant component. 

In order to prove the above claim, let e > be sufficiently small so that (1 — e)A > Ac. We 
use the thinning property to split IIj into two Poisson point processes, 11^ and 11'/, with intensities 
(1 — e)A and eA respectively. Let G[ and G^+i be the RGGs induced by 11'^ and '^[j^i respectively. 
Then with probability 1 — e~ '^'^ ' both G'^ and G^_|_i contain a giant component |28j . We show 
that at least one node from 11/ belongs to both giant components. For any node v of 11/, the 
probability that v belongs to the giant component of G[ is larger than some constant c = c{d). 
Moreover, using the FKG inequality we can show that v belongs to the giant components of both 
G'j^ and G^+i with probability larger than <?. Therefore, using the thinning property again, we can 
show that the nodes from 11/ that belong to the giant components of both G[ and G^,^ form a 
Poisson point process with intensity eAc^, since c does not depend on 11/. Hence, there will be at 

2 

least one such node inside Sn with probability 1 — e~'''^ ", and this stage is concluded by taking the 
union bound over time steps i. 
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Our ultimate goal is to show that if Tperc < t then all nodes of G receive the message being 
broadcast within 2t steps w.h.p. We proceed to the second stage of the proof, and show that the 



tail bound on Tperc from Theorem 1.2 also holds when applied to the finite region Sn defined above. 
Note that all the derivations in the proof of Theorem 1 1 . 2 1 wer e restricted to the cube Q2L, where L 
was defined near the beginning of Section |4j Therefore, it is enough to show that Q2L is contained 
inside Sn (so that the toroidal boundary conditions do not affect the result). But this holds for all 
sufficiently large n since L = 0{t) = 0(log ~^ ' n) while Sn has side-length {n/XY''^. 

The last stage of the proof consists of showing that adding a node u at the origin and calculating 



its percolation time (as we did in Theorem 1.2 ) is equivalent to calculating the percolation time of an 
arbitrary node of G. Note that, by a Chernoff bound, G has at most {l + 6)n nodes with probability 
larger than 1 — e~^^^' for any fixed (5 > 0. These nodes are indistinguishable, so letting p be the 
probability that an arbitrary node has percolation time at least t, we can use the union bound to 
deduce that this applies to at least one node in G with probability at most (1 + 5)np. Let v be an 



arbitrary node. In order to relate p to the result of Theorem 1.2 we can use translation invariance 
and assume that v is at the origin. Then, by the "Palm theory" of Poisson point processes |31j . 
p is equivalent to the tail of the percolation time for a node added at the origin, which is precisely 



d 



Pr [Tperc > t\. Thus finally, using Theorem 1.2 we get p < exp(— ct''+2), which can be made o{l/n) 



by setting C sufficiently large. This completes the proof of Corollary |1.3[ 

Remark: It is easy to see that the above result also holds in the case where the MGG has exactly n 
nodes. The proof above shows that, by setting G large enough, we can ensure Pr [Tbc > i] = o(l/n) 
for the given value of t. Also, it is well known that a Poisson random variable with mean n takes 
the value n with probability p = Q{1/ ^/n). Therefore, for a MGG with exactly n nodes, we have 
Pr[Tbc<t] = 2zi^ = l_o(l/V^). 

6 Some open questions 

We hope our work will promote further mathematical research on mobile networks. Some natural 
open questions related to our results include the following: 

1. What is the tight asymptotic behavior of the tail of the percolation time Tperc? We conjecture 
that the upper bound of Theorem |1.2| can be tightened to match the lower bounds provided 
by Pr [Tdet > i] in Theorem |1.1[ 

2. Let the coverage time Tcov be the time until all points of a finite region S are detected by 
the MGG. How does Tcov behave? 

3. Let Tcomm denote the time until a specific node u is able to send a message to a target 
node V, assuming that the other nodes cooperate maximally to achieve this following the 
broadcast protocol of Section [5] Since Tcomm is bounded above by the time until both u 
and V simultaneously belong to the giant component, following along the lines of the proof 

I I d 

of Theorem 1.2 for the tail of Tperc we can show that Pr [Tcomm ^ i] = exp(— 0(i''+2)). Can 
one show a substantially better upper bound on the tail of Tcomm? 

4. What is the broadest (realistic) class of mobility models (i.e., beyond the Brownian motion 
we consider in this paper) for which our results still apply? 



In our result on the broadcast time (Corollary 1.3), we assumed that messages can travel 
instantaneously throughout each connected component, which is reasonable in many applica- 
tions (where, e.g., transmission through the air is effectively instantaneous in comparison to 
the motion of the nodes). How is this result affected by the assumption that messages may 
only travel a limited number of hops at each time step? 
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A Standard large deviation results 

We use the following standard Chernoff bounds and large deviation results. 

Lemma A.l (Chernoff bound for Poisson). Let P be a Poisson random variable with mean A. 
Then, for any < e < 1, 

Pr [P > (1 + e)A] < exp ( - — (1 - e/3) 
18 



and 

Pr [P < (1 - e)A] < exp ( — — 

Lemma A. 2 (Large deviation for normal). Let N be a Normal random variable with mean and 
variance o"^. Then, for any x > 0, 

a f x^ 



^■■i~^^"^vi;fi'='''v 2.= 
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